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Abstract 

Understanding the superconducting properties of MgB2 is based strongly on knowledge of its electronic structure. In 
this paper we review experimental measurements of the Fermi surface parameters of pure and Al-doped MgB2 using 
the de Haas- van Alphen (dHvA) effect. In general, the measurements are in excellent agreement with the theoretical 
predictions of the electronic structure, including the strength of the electron-phonon coupling on each Fermi surface 
sheet. For the Al doped samples, we are able to measure how the band structure changes with doping and again these 
are in excellent agreement with calculations based on the virtual crystal approximation. We also review work on the 
dHvA effect in the superconducting state. 
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1. Introduction 

The physical properties of MgB2 are strongly 
linked to its unusual electronic structure. This has 
been extensively investigated theoretically and has 
led to a detailed microscopic understanding of the 
physics of this material. Experimentally, there are 
many probes which are sensitive to various aspects 
of the anisotropic electronic structure, although 
there are relatively few which can determine the 
most fundamental parameters accurately. For ex- 
ample, measurements of the resistivity can be used 
to estimate the average mean-free-paths I along the 
two principal crystallographic directions (a and c), 
but cannot provide fc-resolved information. 

Knowledge of these fe-dependent parameters and 
how they vary with doping is important for under- 
standing the superconducting properties of MgB 2 , 
such as T c and the upper critical field H C 2 . In partic- 
ular, a knowledge of how a certain dopant changes 



the band filling and the scattering rates on the var- 
ious Fermi surface sheets can help in optimizing the 
superconducting properties for applications. For ex- 
ample, in order to raise H C 2 it is necessary to lower 
the mean- free-path, especially on the a sheets of the 
Fermi surface, which have a large superconducting 
gap, without depressing T c too much by mixing in 
states from other parts of the Fermi surface where 
the gap is smaller. 

The de Haas- van Alphen (dHvA) effect, is a 
powerful method for determining the Fermi surface 
parameters of metals, providing information about 
the size and shape of the Fermi surface sheets, 
and the scattering rates and quasiparticle effective 
masses on these sheets pp. The data are usually 
interpreted using Lifshitz-Kosevich (LK) theory 
PP. The main experimental constraint comes from 
the fact that dHvA oscillations are only observable 
if the quasiparticles complete a significant frac- 
tion of a cyclotron orbit before being scattered. In 
practice, high magnetic fields and reasonably pure 
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Fig. 1. (color) Calculated Fermi surface of MgB2, with pos- 
sible dHvA extremal orbits (for frequencies <10 kT) indi- 
cated. Panels (a)-(d) show the cr-light hole, cr-heavy hole, 
7r-hole and 7r-electron bands respectively. 

samples with long mean- free-paths are needed. This 
limits the study of doping effects on the electronic 
structure, since doping is usually achieved by par- 
tial atomic substitution which inevitably increases 
scattering and therefore decreases £. The ampli- 
tude of the dHvA signal is reduced from its value 
in a perfectly clean sample by the Dingle factor 
PQ, Rd = exp(-^) = exp(-^gf ). The damping 
is increased with shorter t and increased average 
Fermi wavevector (kp) of the orbit (or equivalently 
as the product of the cyclotron frequency lu c and 
the scattering time r decreases). This can be offset, 
to some extent, by making measurements at very 
high magnetic fields B, although in practice, we are 
still limited to dopant levels of a few atomic percent 
where the increase in scattering is not too severe. 

In this paper we will review studies of the dHvA 
effect in both pure and aluminum doped MgB2, and 
information obtained about the electronic structure 
and superconducting properties. 

2. Electronic structure 

The interpretation of dHvA data is considerably 
aided by detailed comparison with band structure 
calculations. The first calculations of the band struc- 
ture of MgB2 were reported long before the discovery 
of superconductivity in MgB 2 (see Ref. [5]). Shortly 



after the discovery of superconductivity, Kortus et 
al. [3] reported detailed calculations of the Fermi 
surface which specifically address the issue of the ori- 
gin of the superconductivity. Since then there have 
been a large number of other reports which have con- 
firmed and added to this picture (see for example, 
Refs. [4I5I6I7I8Q . 

Generally, the electronic orbits which give rise to 
dHvA oscillations are located at the points where 
there is a local minimum or maximum in the Fermi 
surface cross-sectional area along the direction of 
the applied magnetic field. The observed dHvA fre- 
quency is proportional to this 'extremal area'. In 
most cases these occur at obvious symmetry points 
in the reciprocal unit cell, although the origin of 
some orbits is more subtle. Several papers |7|9|8j 
have reported calculations of the dHvA extremal 
cross-sections and effective masses, and generally 
these are in good agreement with each other. 

In Fig. m we show the Fermi surface of MgB2, 
calculated using the Wien2K package fTO] (see 
later). We have indicated the locations of the pre- 
dicted dHvA orbits [7I9I8| . These have been labelled 
F\ . . . Fq in accordance with previous publications 
[11112113114] . 

3. Experimental details 

The Landau level quantization of the energy lev- 
els of a metal in a strong magnetic field leads to 
oscillations of the density of states as a function of 
the magnetic field and thus to corresponding oscil- 
lations in the physical properties PQ . The most pop- 
ular technique is to measure the oscillations in the 
differential magnetic susceptibility of the sample. In 
practice, the sensitivity of this technique is limited 
by the volume of the sample. In the case of MgB2, 
the available highest quality single crystals are rel- 
atively small (<0.05mm 3 ), so better signal to noise 
can be obtained by measuring the oscillations in 
the torque, using a piezo-resistive cantilever of the 
type developed for atomic force microscopy [15j . All 
the results we will show here were obtained by this 
method. The details of the method are described in 
Ref. [12]. 

4. Results for pure MgB2 

The first report of dHvA oscillations in MgB2 was 
by Yelland et al. [11) . In this report, which used mag- 
netic fields up to 18 T, three dHvA frequencies (Fi, 
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i<2, and F3) were observed. F\ and F2 correspond to 
the extrema of the smaller (light hole) a tube, and 
F3 is located on the electron-like ir sheet. In a later 
study on crystals produced from isotopically pure 
10 B, Carrington et aL[13] observed oscillations from 
all four sheets of the Fermi surface (F\ . . . Fq in Fig. 
[T]), using fields up to 33 T at the Tallahassee high 
magnetic field laboratory. 

The observed dHvA frequencies as a function of 
field angle are shown in Fig. [2l Results for two dif- 
ferent samples are shown. Sample B was grown from 
natural (mixed isotope) boron whereas sample K 
was the pure 10 B sample. In crystal K signals from 
all four sheets of the Fermi surface were observed, 
however, in crystal B no oscillations arising from the 
larger (heavy hole) a sheet were observed even in 
these higher fields. This is almost certainly because 
of the longer mean- free-paths on the a sheets of crys- 
tal K relative to those in B (see Table [l| . Note the 
opposite trend of £ on the ir sheets of the two crys- 
tals. In general, we find that the mean-free-paths on 
the a sheet are not well correlated with those on the 
7r sheets. This will be discussed later. 

The dHvA frequencies, extrapolated by polyno- 
mial fits to the symmetry points F°, are given in 
Table [TJ We find very good reproducibility between 
samples from different sources, and the dHvA fre- 
quencies agree to within 30 T or 0.06% of the basal 
area of the first Brillouin zone (^^j = 50.2 kT). 
This is evidence against there being any significant 
Mg deficiency in the samples. In general, the agree- 
ment between the measured F values and those 
predicted by theory is very good although not per- 
fect. To investigate the size of the discrepancy we 
apply a rigid shift to the calculated energies of the a 
and 7r bands. Since the dHvA band mass is defined 
as ms = l^fjj, the necessary band shifts are given 
by AE = ^AF (where AF = F LBA - F exp ). For 
the a sheet orbits (1,2,5,6) the average shift is 83±4 
meV, whereas for the n sheet orbits (3,4) it is — 61±5 
meV. The volumes of the c-axis tubes are propor- 
tional to the average of the two extremal areas [8] 
and these are both ~ 16% smaller than the calcu- 
lations |I6j . implying a corresponding difference in 
the number of holes in these two tubes [8] . 

In Fig. [3]we compare the frequencies observed ex- 
perimentally with the predictions based on the LDA 
band structure. Here we have re-calculated the band 
structure using the Wien2K package [TOj (using pa- 
rameters similar to those of Kortus et al. [3] ) on a 
very dense fc-mesh (95000 points in the full Brillouin 
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Fig. 2. Observed frequencies versus field angle as the samples 
were rotated from H\\c to (approximately) H\\a. The solid 
lines are polynomial fits , and the dotted lines are guides to 
the eye. 

zone) and extracted the dHvA parameters using an 
automated numerical search algorithm. The bands 
have been rigidly shifted to give the best agreement 
at the symmetry points as described above. The 
agreement is excellent, showing that the LDA cal- 
culations (with an applied small rigid shift) give the 
correct Fermi surface. 

In MgB2 the dominant source of mass enhance- 
ment is the electron-phonon interaction |5I6I9| . If 
we assume this is the only source of enhancement 
we can estimate the electron-phonon coupling con- 
stants A cp , from the measured quasiparticlc effective 
masses m*, A cp = m* /ms — 1 • The results (Table[T]) 
show that the values of A ep on both the a sheets are 
approximately a factor three larger than those on 
the 7r sheets, in agreement with theoretical predic- 
tions. This is a key factor in the multi-band nature 
of superconductivity in MgB2, and results in a much 
larger superconducting gap on the a sheets than on 
the 7r sheets. 

The experimental values of A cp are generally 
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Table 1 

Summary of dHvA parameters for both samples (K and B) along with the theoretical predictions (Th) [9] 



Orbit 



Crystal K 



[T] [T] [T] [meV] [A] 



m ba : 



. cp 
A cxp 



A Th 



Crystal B 



F° 

exp 

[Tl 



I 

[A] 



A cxp 



551 730 +179 +83 550 0.548 + 0.02 0.251 1.18 + 0.1 1.25 

1534 1756 +222 +82 900 0.610 + 0.01 0.312 0.96 + 0.03 1.25 

2705 2889 +184 -67 570 0.439+0.01 0.315 0.40 ± 0.03 0.47 

576 458 -118 -56 - 0.31 + 0.05 0.246 0.31 + 0.1 0.43 

2971 3393 +422 +79 390 1.18 + 0.04 0.618 0.91 + 0.07 1.16 

1180 1589 +409 +87 - 1.2 + 0.1 0.543 1.2 + 0.2 1.16 



546 380 0.553 + 0.01 1.20 + 0.04 

1533 580 0.648 ± 0.01 1.08 ± 0.03 

2685 680 0.441 ± 0.01 0.40 ± 0.03 

553 - 0.35 + 0.02 0.42 + 0.08 
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Fig. 3. (color) Comparison of the observed dHvA frequencies 
with extremal FS areas extracted from the calculated band 
structure (A=crystal K, »=crystal B). Calculated frequencies 
are for an in-plane angle <f> = 12°, to match crystal B. 

slightly smaller than the theoretical ones (even 
when the a and it band shifts, mentioned above, are 
taken into account p2] ■ One possible reason for this 
is that these theoretical values do not take into ac- 
count phonon anharmonicity, which has been shown 
to reduce the average value of A by around 20% [T7] . 

In sample B, two additional dHvA frequencies 
were observed which we did not see in crystal K 
(these are labelled F 3 * and F 3s ). When the field is 
aligned in-plane, close to the a-axis [100], the orbit 
-F3 can be seen, but in addition there are also two or- 
bits in the symmetrically equivalent section of Fermi 
surface running along the [110] and [010] directions, 
inclined at 60° to the main F3 orbit. The origin of 
Fq s has not yet been conclusively identified. 

As far as we know, there has been only one other 



report of dHvA measurements in MgB2 by another 
group. Isshiki et al. also used a piezoresistive can- 
tilever method to measure crystals grown by an en- 
capsulation technique [14] . Their results are very 
similar to those shown above for crystal B (Fig. \2§- 
They observe all the same frequencies (including F 3 * 
and F 3s in Fig. [2]), but do not see any signals from 
the heavy hole a sheet (F 5 and F 6 ), most likely be- 
cause of a lower mean-free-path in these samples. 
The frequencies at the symmetry points {F^) are 
within 20T of those in Table Q] except F4 which is 
46 T lower ITS]. 



5. Results for Al-doped MgB 2 

Studies of the effect of atomic substitutions in 
MgB2 are important for applications and may 
provide a route to improve its superconducting 
properties. It also allows us to test theoretical 
understanding of the material. The two elements 
which substitute most readily in MgB2 are Al and 
C which replace Mg and B respectively, to give 
Mgi_ a; Al a: (Bi_ y Cj / )2. Both dopants add electrons 
to the material, and cause T c to decrease in a simi- 
lar way [19]. However, the effect of each dopant on 
H C 2 and resistivity is very different. Al doping gives 
a moderate increase in the residual resistivity and 
causes only a small change in H C 2 and its anisotropy. 
On the other hand, carbon doping dramatically in- 
creases both the resistivity and H C 2 [20121) . This 
is understandable as the a band has most weight 
on the boron plane, whereas the tt band has weight 
on both the boron and magnesium planes. It is 
therefore to be expected that replacing boron with 
carbon should strongly increase the scattering rates 
on both the a and 7r bands, but replacing Mg with 
Al would mainly affect the 7r band. 
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Table 2 

Calculated (LDA) 9 and measured dHvA frequencies [F° = 
F(6 = 0)] for pure [13] and Al doped MgB 2 . AE is the rigid 
band shift needed to bring the theoretical values in line with 
experiment. 
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Fig. 4. Observed dHvA frequencies versus angle for the two 
Al doped samples AN215 (■,♦) and AN217 (0)- Dashed 
lines are F(0) for pure MgB2, solid lines are pure F(0) data 
scaled to fit the Al doped data. 



dHvA studies can test theoretical predictions by 
measuring sheet-specific scattering rates and by de- 
termining quantitatively how doping changes the 
volumes (or cross-sections) of Fermi surface sheets 
and the effective masses of quasiparticles. The latter 
can indicate the doping-dependence of the electron- 
phonon coupling constants. 

Attempts were made to detect dHvA oscillations 
in both Al and C doped samples with similarly re- 
duced T c [55]. The C doped samples had T c values 
of 35.7 K and 34.5 K, and nominal C contents of y = 
3% and y = 4% respectively. Their residual resis- 
tivities p(T c ) were ll/zf2cm and 30 ^.Q cm, substan- 
tially larger than for pure MgB2 which typically has 
p(T c ) ~ lfiCl cm. No dHvA signals from either the a 
or 7r sheets were seen in either C doped sample. The 
two Al doped samples (AN215 and AN217) both had 
T c ~ 33.6 K (as determined by heat capacity mea- 
surements), and Al contents of x = 7.9 ± 0.4% and 
x = 7.4 ± 0.4% (as determined from their c-axis lat- 
tice parameters [53]). The residual resistivity of the 
Al doped samples was found to be ~ 4.5 /Ltficm. 

In Fig. [5] we show the observed dHvA frequen- 
cies versus field angle, 6. Only signals arising from 
the small (light hole) a sheet were observed. This 
sheet has the smallest cross-section (or Uf) so the 
increased scattering (or decreased €) has the small- 
est effect on the Dingle damping term. The values 
of Ff and F® are the same in both Al doped crys- 
tals, and are 10 — 20% smaller than the correspond- 
ing values for pure MgB2 (see Table [5]) . Studies on 
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many different samples of pure MgB2 (from differ- 
ent sources) have revealed remarkably reproducible 
dHvA frequencies and so the measured reduction in 
cross-section is significant. Within a rigid band shift 
model our results imply that at this doping the a 
bands have shifted to lower energy by ~70 meV rel- 
ative to the experimental results for the pure sam- 
ples. From the calculated volume of these sheets we 
find that there are 16 ± 2% fewer holes in 7% Al 
doped samples compared to pure MgB2 . 

This simple analysis was improved considerably 
by J. Kortus [19)22] who calculated the effect of Al 
doping on the band structure using the virtual crys- 
tal approximation ( VCA) . Here the effect of doping 
with Al is simulated by replacing the Mg atom with 
a virtual atom with charge Z = xZal — (1 — s)^Mg- 
These calculations [22] give the following depen- 
dence of dHvA frequencies on Al content x, F\ = 
540 - 1820a;, F 2 = 1530 - 2050x (here the bands 
have been rigidly shifted to give agreement with the 
experimental data at x = 0, i.e. for pure MgB2). 
The observed frequencies (-Fi and F%) correspond to 
a doping of a; = 7.5 ± 1% and 8.4 ± 1% respectively. 
These values compare favorably to the x values de- 
duced from the c-axis lattice constant of the same 
crystal (x — 7.7 ± 0.4%), and so we conclude that 
the VCA results accurately describe the band filling 
effect of Al doping in MgB 2 . 

The VCA calculations are also able to predict 
the decrease in T c expected from the doping. For 
each doping the density of states, phonon frequen- 
cies and electron-phonon coupling parameters were 
calculated. Close to x = the calculations fH] pre- 
dict dT c /dx = -0.50 K/%, and hence for x = 7.7% 
we expect a T c reduction of 4.0 ± 0.3 K. The actual 
T c reduction in our samples was 3.9 K, in excellent 
agreement with the calculations. Note that these cal- 
culations do not include the effect of the increased 
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scattering, so it would appear that, at least for this 
doping level, the scattering has a minimal effect on 
T c . 

A further test is to look at the decrease in the mea- 
sured electron-phonon coupling strength with dop- 
ing. For the Al doped samples, the only orbit where 
we have sufficient accuracy to make a meaningful 
comparison is F%. For this orbit we find that A ep is 
10 ± 5% smaller than for the pure case (the quoted 
error includes the measured variation in A ep between 
different pure samples) 22J. In calculating this we 
have taken into account the small change in the band 
mass with doping. Theoretically, it is found that A ep 
does not vary substantially within each a sheet and 
only varies by ~ 7% between the two a sheets [9] , so 
the observed reduction in A ep measured on orbit F% 
is expected to be representative of the changes on 
the other a sheet. For x — 7.7%, the VCA calcula- 
tion [19122] gives A£ p = 1.16, which is - 5% smaller 
than for x = 0. This is consistent with our observa- 
tions. 

The mean-free-path, estimated from the Dingle 
term i?£>, was I = 270 A for orbit F%. This is ap- 
proximately half of that found for pure samples of 
MgB2 (grown with natural mixed- isotope boron). 
The signal from the electron-like 7r sheet (frequency 
F-i) is normally large for fields oriented near the 
plane, [11113] and as we have not observed this or- 
bit in either Al doped sample we conclude that the 
mean-free-path on this sheet must be reduced by a 
somewhat larger factor (> 3). 

Using this information we can make some sim- 
ple estimates of the effect of Al doping on H C 2- In 
the clean limit, fJ,oH C 2 = $o/(27r£ 2 ), where £ can 
be estimated from the BCS coherence length £o = 
Hvf /(""A). At low temperature the a sheets domi- 
nate p4 25 26 27] and the anisotropy of H c2 (jh c2 ) is 
determined by the anisotropy of vf on the a sheets, 
i.e., lHc2 = H%/H% = \/(v F ,a)/(vF,c)- In the dirty 
limit, the same result is found [25] provided the scat- 
tering is isotropic. LDA calculations show that for 
x < 0.3, [26128] 

y/(v Fta )/(v FtC )^5.9-7.5x. (1) 

Experimentally, we find for pure MgB2, ^oH^ = 
3.5±0.2Tand7// c2 = 5. 5±0.2, at T ~ 0.3K, whereas 
for our Al doped samples (x = 7%), HqH^ = 3.5 ± 
0.1T and 7# c2 = 4.6 ± 0.2. These values are in 
line with other studies [29126127130121131] . Eq. Q} 
slightly overestimates the anisotropy; however, the 



agreement improves if we take into account the fact 
that pure MgB2 has smaller a tubes compared to the 
LDA calculations. The hole doping of the a sheets 
for pure MgB 2 is roughly equivalent to that in the 
LDA calculations with x ~ 11%, so the predicted 
Jh C 2 nas a value of 5.2 for the pure case and 4.6 for 
the Al doped (x=7%) case. These latter values are in 
good agreement with experiment. In addition, using 
the values of A and vf for the a sheets, £o — 130A. 
This is less than I on these sheets for both the pure 
and x=7% samples and hence the lack of significant 
change in (j,qH C 2(T = 0) is also understandable. In 
other words the a bands for both the 0% and 7% 
Al-doped samples are moderately clean (£q < t) . 

6. Angle dependence of dHvA amplitude 

The angular dependence of the dHvA ampli- 
tude is often used to extract information about 
the strength of electron-electron interactions. The 
Lifshitz-Kosevich expression describing dHvA os- 
cillations in a paramagnetic metal contains a factor 
Rs = cas[^irg(l + S)m-Q / m e \ which arises from the 
phase difference between dHvA oscillations from the 
spin-split Fermi surfaces pQ . Here, S is the Stoner en- 
hancement due to the electron-electron interactions 
that enhance the Pauli spin susceptibility. At certain 
field angles, oscillations from the spin-split Fermi 
surface sheets may destructively interfere; Rs then 
vanishes and a 'spin-zero' is observed in the dHvA 
amplitude. Measurement of the field-angle 6 SZ for a 
spin- zero can then be used to determine the Stoner 
enhancement S — (2/g)(n+ ^)m e /m^,{9 sz ) — 1, as 
long as the band-mass 77Jb(#), n and g are known. 
n can often be deduced from the frequency of the 
'spin-zeros' as a function of angle; tob(^) is taken 
from LDA band structure calculations. 

Fig. [5] shows the dHvA amplitude versus field an- 
gle for orbits F\ , F2 and F3 . The dips in the ampli- 
tudes of Fx and F 2 are due to spin-zeros [T3] ; fits to 
the LK expression including i?g (solid line) are in 
excellent agreement with the data. The Stoner en- 
hancements are estimated to be 5(Fi) = 0.10 and 
"5(^2 ) = 0.14 [32]. These are approximately a fac- 
tor 2 lower than the values from band structure cal- 
culations [9]. The dHvA values are more consistent 
with conduction electron spin resonance measure- 
ments of the spin susceptibility [33] , which imply no 
enhancement (S = 0) to within the experimental 
error of ±15%. 

A sharp minimum in the amplitude of the F 3 orbit 
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Fig. 5. Top panel: dHvA torque amplitude versus angle for 
frequencies Fi, Fi and F3 (17.8 < B < 18.0 T). For frequen- 
cies i*i and F2 the solid lines are fits to the LK expression 
including the spin-splitting factor i?g . For F3 the dotted line 
is a guide to the eye. Bottom panel: F(k^) = hA(k^) /2ne 
for certain field angles, from LDA calculations. For fields 
close to the boron plane, both minimal and maximal cross- 
sections are present, but for 9 < 65° there is only one ex- 
tremal orbit. Inset: the amplitude of F3 fitted using a two 
frequency model (solid line) and calculated from the LDA 
band-structure (dashed line, see text). 



on the 7r band sheet was also observed (at 9 « 74°, 
4> = 0°). This was initially interpreted as a spin-zero 
[11112113] . although this required an anomalously 
large Stoner enhancement. A later detailed dHvA 
study [33] revealed that the dip is in fact due to a 
small warping of the nearly cylindrical part of this 
sheet which can be seen in LDA calculations (see 
Fig. [5]) . This warping means that there are two ex- 
tremal dHvA orbits which have very similar frequen- 
cies for 9 > 74° . For lower 9 only a single frequency 
is observed. Fits of the observed angular dependence 
using a model with two discrete frequencies sepa- 
rated by AF ~ 1.4(0 — 70°) T can reproduce a sharp 
minimum (see Fig.[S]), but are not completely satis- 



factory [34] . It is well known that the extremal or- 
bit approximation used in the LK formula breaks 
down when there are significant regions of the Fermi 
surface over which F(kn) varies by < B, so that 
changes in the dHvA phase cannot be treated as 3> 
2t:. The effect is predominant in metals with a quasi- 
two-dimensional Fermi surface |35|36j . but it is also 
significant in the nearly cylindrical sections of the 
electron-like it band Fermi surface sheet in MgB2. 

Fig. ED shows F(k\\) = 7L4(fe||)/27re [where A(k\\) 
is the Fermi surface cross-sectional area] , calculated 
from LDA results, at several field angles. For angles 
> 65°, F(k\\ ) has two distinct extremal values. The 
contribution of this finite Fermi surface region to 
the dHvA amplitude can be calculated by evaluating 



2rrF(k ll 
B 



over a suit- 



the integral 1(B) = J dk\\ sin 

able range of kn ■ The inset of Fig.[5]shows the dHvA 
torque amplitude calculated in this way, taking ac- 
count of the angle dependence of Rd, Rs an d &F/&6. 
The calculated amplitude has a series of sharp min- 
ima which resemble the single minimum seen in the 
experimental data. The reason for the discrepancy 
in the number and angular location of the minima 
is most likely due to a small inaccuracy of the LDA 
calculation: a warping approximately 30 T smaller 
than that predicted by the LDA (at 9 = 90°) would 
place the first minimum at the experimentally ob- 
served angle (other minima would be shifted to an- 
gles where strong damping would prevent their ob- 
servation). This level of error in the LDA calculation 
is entirely plausible because the value predicted for 
F 3 is ~ 200 T too high before applying a rigid energy 
shift. The absence of a true spin-zero for 9 > 50° im- 
plies a Stoner enhancement S < 0.22, in agreement 
with other estimates. 

7. Damping in the superconducting state 

In strongly type II, high H C 2, superconductors, 
quantum oscillations can often be observed in the 
superconducting state, i.e., at fields where supercon- 
ducting vortices penetrate the sample (see for ex- 
ample, Ref. [37]). In pure MgB 2 , [i H c2 - 18 T for 
fields oriented close to the boron planes, and for field 
angles 9 > 65° we are able to observe dHvA oscilla- 
tions from the ir band orbit F3 below H C 2 [38]. Fig. [6] 
shows an example of torque data obtained close to 
H C 2- Quantum oscillations are superimposed on the 
background torque due to superconductivity, which 
arises from the anisotropy of the coherence length 
and pinning [39] . As is generally found in other mate- 
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Fig. 6. Torque measured near T c with 9 = 69.7° and ~ 5°. 
At this angle H c 2 ~ 11.5 T, as estimated from the rest of the 
torque loop. The substantial irreversibility above this field is 
probably due to surface superconductivity. Note that this is 
part of a minor hysteresis loop; in large samples the torque 
arising from superconductivity is too large to measure over 
a complete loop. 

rials, the dHvA frequency remains identical to that 
in the normal state but the oscillations suffer an ad- 
ditional attenuation beyond the factors accounted 
for in the conventional Lifshitz-Kosevich treatment. 

The survival of quantum oscillations deep into the 
superconducting state is surprising for several rea- 
sons: (i) the vortex state is characterized by an in- 
trinsic field inhomogeneity on the scale of the Lon- 
don penetration depth, (ii) pinning of vortices leads 
to a vortex density gradient and thus field inhomo- 
geneity on the scale of the sample, and (iii) the open- 
ing of the superconducting gap modifies the low en- 
ergy quasiparticle spectrum, reducing the number 
of carriers available to undergo cyclotron motion. 
Despite this, oscillations can sometimes be seen at 
fields as low as 0.2 H C 2 [40] . A great deal of theoret- 
ical work [41] (see Ref. [42] for a review) has been 
performed to establish whether the origin of the os- 
cillations is the same as that in the normal state, and 
to understand the mechanisms which attenuate the 
oscillations. Many models attribute the attenuation 
to the reduction in the (spatially averaged) density 
of states at the Fermi energy in the vortex state [43] . 
which is related to the orbit-averaged superconduct- 
ing gap. This is particularly interesting as the size 
and fc-dependence of the additional damping might 
be used to determine the fc-dependence of the super- 
conducting gap. Rigorous experimental tests of the 
various theoretical models for damping of quantum 
oscillations in the superconducting state have been 
hard to perform because such observations are gen- 
erally made in exotic materials where even the nor- 



mal state is not well understood: in some cases only 
a fraction of Fermi surface is detected (even in the 
normal state) and the nature of the effective mass 
renormalization is unknown. Our knowledge of the 
electronic structure and superconducting properties 
of MgB2, make this an ideal material for performing 
such a study. 

A common way of isolating the additional atten- 
uation is by dividing out other field dependent fac- 
tors from the original LK expression (see Ref. [38] 
for details). This additional damping factor, R sc , is 
shown in Fig.[7]for MgB2. The results were found to 
be similar at all measured angles, apart from chang- 
ing the field values at which the attenuation begins. 
The attenuation has a rather rounded onset, but is 
quite severe, corresponding to an order of magni- 
tude, at 0.9 H c2 . 

The value at which attenuation becomes signifi- 
cant is substantially lower than the point at which 
the torque becomes irreversible. In fact, as shown in 
Fig. irreversibility extends to much larger fields 
than is expected from the rest of the superconduct- 
ing torque loop. This is most likely due to surface su- 
perconductivity [14144] , although some calculations 
suggest that fluctuation effects may be important in 
this region [45] . 

The field dependence of R sc was found to be in- 
dependent of sweep direction, i.e., unaffected by the 
irreversibility of the background torque, and was re- 
produced in other samples with different degrees of 
vortex pinning. This suggests that the signal is unaf- 
fected by inhomogeneity of the the vortex distribu- 
tion due to pinning. Following the arguments in Ref. 
[3"7] , it is clear that the small scale field inhomogene- 
ity due to a homogenous vortex lattice also does not 
produce significant attenuation of the dHvA signal 
in MgB 2 . 

In another study of this effect in MgB2 [14] , an ac 
field was superimposed on the main field and modu- 
lations in torque were detected. If vortices are pinned 
at the sample surfaces, the ac field modulation does 
not penetrate the bulk of the sample and an addi- 
tional damping of the observed signal results. The 
dc field torque method used in the work described 
here [38] is not affected by this problem. 

The damping of the dHvA signal is closely linked 
with the onset of superconductivity in the bulk of 
the sample. In the present case, additional attenu- 
ation seen in the dHvA signal is attributed to the 
finite gap present on the ir band, which is present 
even in fields close to H C 2- Very few probes of the 
superconducting gap exist in such high fields, and 
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Fig. 7. Additional attenuation of the oscillatory torque signal 
in MgB2 in the superconducting (sc) state. For 9 = 69.7°, 
data taken on up sweeps (o) and down sweeps (•) are shown. 
Data are shown for in-plane rotations <j> = 5° and <j> = 24°. 
Due to fine details of the 7r-band Fermi surface sheet (see 
text), measurements at <f> = 24° allow R BC to be determined 
over a wider field range. 

in fact it was previously suggested that the super- 
conductivity in the n bands was quenched at much 
lower fields [46] . 

In Fig. two models are applied to the measured 
damping. The zero field superconducting gap A 
and H C 2 are free parameters. An additional param- 
eter is included to account for the broadening of the 
transition [47]. The first model [48] is based on the 
calculation of the average density of states in the 
vortex lattice 43J. This requires A ~200-300 K 
in order to describe the data. An alternative model 
[4T)] , simply based on the zero field (and hence fully 
gapped) density of states, requires a value of Ao ~ 
30-60 K, closer to that expected in MgB 2 . Compared 
to results in other materials [37] this level of agree- 
ment is not unusual. While this may indicate that 
this effect has yet to be fully understood, we also 
note that there are some limitations in this analysis. 
Firstly, data are only available at high fields, requir- 
ing a substantial extrapolation of the field depen- 
dence of the gap back to the zero field value. This 
field dependence is non-trivial in the multi-band case 
[50] . In addition, the meaning of the gap parame- 
ter in each model can be different, i.e., it represents 
an effective gap, which may be quantitatively dif- 
ferent from the zero field value. In the light of this, 
it would be interesting to observe dHvA signals be- 
low H C 2 on both a and 7r bands in MgB2 to confirm 
that the relevant parameter is the sheet-dependent 
(orbit-averaged) gap. This would require the dHvA 
signal to survive down to fields of order HqH^ — 



3.5T. We estimate that this would require samples 
with a sheet mean-free-paths approximately three 
times larger than we measured to date [38] . 

8. Conclusions 

In conclusion, we have seen that in general, our 
measurements are in excellent agreement with the 
predictions of the LDA band structure, including 
the strength of the electron-phonon coupling on each 
Fermi surface sheet. The minor discrepancies which 
remain can be explained with small rigid band shifts. 
Our work on Al substitution confirms that adding 
Al dopes electrons into the structure, and reduces 
the size of the a sheets. The size of the reduction and 
the size of the T c decrease is quantitatively explained 
by the VCA calculations. 

Although our understanding of the dHvA data 
obtained until now is rather good, the data for Al 
doped samples show that there is scope for further 
studies of crystals lightly doped with other elements. 
These would give important insight as to how the 
mean-free-paths on various orbits and the strengths 
of the electron-phonon coupling can be altered in a 
controlled manner with alloying and give deeper un- 
derstanding as to how T c and H C 2 can be optimized 
for applications. 

The dHvA effect in the superconducting state, 
however, is still not fully understood. The existence 
of two distinct gaps in MgB2 gives us an important 
new parameter with which to study this effect - espe- 
cially if in the future higher purity samples become 
available so that dHvA oscillations can be measured 
at fields less than /io-ffj,'^ — 3.5T and the supercon- 
ducting gap on the a sheet can also be probed. 
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